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Abstract. We revisit the directionally optimal data streams of LISA first introduced in Nayak etal. 
It was shown that by using appropriate choice of Time delay interferometric (TDI) combinations, a 
monochromatic fixed source in the barycentric frame can be optimally tracked in the LISA frame. 
In this work, we study the beaming properties of these optimal streams. We show that all the three 
streams v+,x,o with maximum, minimum and zero directional SNR respectively are highly beamed. 
We study in detail the frequency dependence of the beaming. 
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1. INTRODUCTION 

The space-based gravitational wave (GW) mission LISA [1] — Laser Interferometric 
Space Antenna — consists of three identical space-crafts forming an equilateral triangle 
of side 5 X 10^ km following heliocentric orbit trailing the Earth by 20°. The plane of 
LISA makes an angle of 60° with the plane of the ecliptic. Each space-craft completes 
one orbit around the sun as well as in LISA plane in one year. The mission is aimed 
at detecting and analyzing the low frequency GW signals in the frequency range of 
0. ImHz — IHz. The astrophysical sources for the LISA include galactic binaries, super- 
massive black-holes (BH), extreme mass ratio inspirals, intermediate mass BHs. 

Due to LISA'S rotational as well as orbital motion, a fixed source in the barycentric 
frame appears to follow a specific track in its sky. This introduces amplitude modu- 
lation, frequency modulation, and phase modulation in the 6 Doppler data streams. A 
large number of interferometric configurations, having different frequency and angular 
response, can be constructed from these data streams which makes LISA not just a single 
detector but a network of interferometers. 

The choice of the combination of the data streams depend on the which question one 
wishes to address. In Ref. [2], we addressed the question of directional optimality in 
LISA i.e. for a given sky location, which LISA data stream gives maximum SNR. It 
was shown that the directional optimality condition gives three data streams: (1) v-|- - 
with maximum directional SNR, (2) Vx - with minimum directional SNR, (3) vq - with 
zero directional SNR. It was shown that the data stream v-|- optimally tracks the source 
motion (in the LISA sky) of a fixed source in the barycentric frame. Here, tracking 
involves appropriate choice of data combinations (switch combinations as source moves) 
which gives maximum SNR in that direction. Tracking known monochromatic binaries 
with such streams could give information about the source distance, polarizations. For 
an unknown distant source, it would amount to 'looking' in a specific direction in a 
particular frequency band. 



In this work, we study the beaming property of the optimal data streams (in particular, 
V-I-) for a monochromatic source tracked for a year. We show that Vj^-.v^andvo are 
beamed, i.e. they are sensitive towards the tracking direction. The beam-width depends 
on the frequency under consideration. We study the nature of this dependence. 

The paper is organized as follows: In Sec 2, we review the TDI. In the first half of Sec. 
3, we summarize the main results of Ref.[2] pertaining to the directionally optimal TDI 
streams. In the later half of Sec 3, we discuss the beaming properties of the directional 
streams followed by conclusion in Sec. 4. 



The 6 LISA Doppler data streams where m = 1,2,3,(7 = ±, [see Fig.l^] are 
obtained by letting the laser beams from each space-craft to travel towards two other 
space-crafts and are beaten with the on board laser, m corresponds to the arm index and 
(— )+ indicates the laser beam traveling in the (anti-)clockwise direction. 



These Doppler data streams contain the phase fluctuation noise of the master laser 
(Av/vo — 10^^^ / \/Hz) which is several orders of magnitude higher than the LISA 
designed sensitivity level of h ~ 10^^^. Using TDI technique, W^^ can be combined by 
appropriately delaying them with time delay operators to construct laser noise-free data 
streams Ref. [3]. The TDI scheme is based on the principle of sagnac interferometry — 
the light-beam is split and made to travel along two paths in opposite direction of equal 
length which are then subtracted to remove the laser noise. In TDI, equal path difference 
is achieved by adding Doppler-streams with appropriate time delay operators which are 
then added (or subtracted) to obtain laser noise free data stream. 

Any laser noise-free data combination can be written as = Lm^i LcT=±PmaW'^, 
where p^a are polynomials of time-delay operators E,n : E,nCj{t) = Cj\t — L^) and Cm 
are the laser phase noise fluctuation at the space-craft m. The Sagnac TDI combinations 
{ a , /3 , 7, ^ } in terms of {p,„+ , pm- } are given by ^ 



a = {£2, l,£i£2,-£3, -£1^3,-1} /3 = {£2^3,^3, 1,-1,-^1, -^1^2}, (1) 

7 = {l,E,E3,Ei,-E2E3,~l,-E2} C = {^3,£i,£2, -^2, -^3, -^i} • (2) 



1 Wl = t/^W| ^U\WI = U^,Wl = -V^,Wl = -V^,Wl = -V'of [3] 



This set forms a generator set for an algebraic module of all laser noise free data combinations Ref. [3]. 



2. TDI DATA STREAMS 




FIGURE 1. The LISA constellation 



The ^ combination is termed as symmetrized sagnac due to its symmetric structure 
and is insensitive to GW at low frequency f « 1/L = 60 mHz (L = 16.7 sec). 
The noise vector for any combination D is given by 

AT^ = {2Vs^{p^+ + il^+)a^{p^-+l^^-),Vs^P^+,Vs^P^-) (3) 

where the polynomials iim± are defined as 113 = {E^p3+ — P3 V2 = — /i2+. The rest 
Hmcy can be obtained by cyclic permutations. The 5^^ = 2.5 x 10^ (// lHz)^^Hz^^ and 
the 3"^"^ = 1.8 X 10^^^(//lHz)^Hz^^ are the one-sided power spectral densities (PSD) 
of the proof-mass noise and optical-path noise respectively [1]. In the frequency domain 
£'m = exp(z^2L„). 

For simplicity, throughout this work, we assume the three arms of LISA to be equal 
i.e. = L. This helps in simplifying the expressions of TDI streams and are exact for 
low frequencies. However for higher frequencies, the above simplification leads to small 
discrepancies. Nevertheless, one can easily extend this for unequal-arm interferometry. 

A set of TDI data streams which diagonalizes the noise covariance matrix A'^^^) • A^^*^^ 
are [4]: y(i) = (a + /3 -27)/v^,y(2) = (j8 - a)/v^,y(3) = (^a + l5 + r)/V3. 7(1), 
and Y^^^ are known as E,A and T in the LISA literature. 

3. DIRECTIONAL DATA STREAMS 

As mentioned in the introduction, the three F^^^'s give different frequency as well 
as angular response. A large number of TDI data streams 0^, (^L)y(^) with 

different angular and frequency responses can be constructed from them. The choice 
of a combination depends on the question one wishes to address. In Ref. [2], we asked 
the following question: If one wants to observe a particular sky location {0s,<^b}, in 
a given frequency bin, which X)ot(/)i'^^^ TDI data stream would be optimal (maximum 
SNR)? We first briefly summarize the results of Ref. [2]. 
The GW response of Y^'^ expressed in frequency domain is 

h^'^ (a) = F^^ {Q.)h+{a) + (a)/ix (n) , (4) 

where, the antenna pattern functions F^^^ in terms of the transfer function of each 
Doppler data stream is 

(5) 

m,a 

where w is the direction vector to the source in the LISA frame, /S.^m = QLm/2, 
= • n appears in the accumulated phase due to GW oscillation as the laser 
travels from one space-craft to another (anti-)clockwise direction, dm is the normal 
vector from LISA centre to each arm m and ^m;+x are the responses of each arm to 
the two GW polarizations. 

The SNR maximization for a particular direction is a constrained optimization prob- 
lem. The data combinations are obtained from the eigen-vectors of the SNR squared 



matrix (averaged over the polarizations) 



P(j) 

where f^!\ - "^^^'^ • 



J+ J + {J) +/x /x(7) 



(6) 



^0'f+,x/'^(/)' ^0 '•^^ average signal energy over the GW polarization 
and is the noise PSD of Y'^^\ The eigen- values are the instantaneous squared SNR 
for the optimal data streams. The 3 eigen vectors are given by: 



C+f++Cxf. 
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4/4 



■C+/x, Vo = /+ x/x 



(7) 



where the coefficient Cx 
1 



/x -Z+and 



I/+I - l/x P + V (I/+P - l/x +4I/+ • /x 



snr I 



l/x 



(8) 



Note that the orthogonal pair of {V+,Vx } is obtained by the linear combination of /+ 
and fx and hence lie in {f+,fx } — polarization — plane. The 2x2 matrix which trans- 
forms {f+,fx } to {V+, Vx } is traceless-hermitian. Naturally, the direction orthogonal to 
this plane contains no signal, i.e. a data stream obtained from Vq (orthogonal to {f+,fx } 
plane) is a null stream. Hence, the triplet {V+^Vx^Vq} gives complete directional in- 
formation of the GW signal. The tracking coefficients for the optimal data streams are 



(/)* 



(/)+ 



y(/) 



gives maximum 



^(/)+,x,o = ^+ X o/"(/)- summary, the data stream (i) v 
directional SNR i. e 
e. snrx and (iii) vq = 

and v+ tracks the + polarization while Vx tracks the x polarization of GW, hence their 
subscripts [2]. 

(a) (b) 



snr+, (ii) Vx = a^j^-^Y^^^ gives the smallest non positive SNR i. 
(^{i)oY^'^ gives zero directional SNR. At low frequency, Cx = 
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FIGURE 2. (a) Source coordinates {Xi{t),YL{t),Zi{t)} in LISA frame for a fixed source at (pB 
200° for various 9b- (b) LISA Sensitivity for various data streams while tracking source at {9b -.^b} 
{60°, 200°} and //o = 1, r = 1 year 



3.1. Source tracking and Integrated SNR 



Due to LISA'S motion, a fixed direction wb — {0b, 0s} in barycentric frame appears 
to follow a track {0^(05 0^(0} LISA'S sky, see Fig. 2(a). We track the direction 



and obtain the integrated SNR as follows: We choose the optimal data combination 
pertaining to the source direction {0L(r ; ws), wb)}, at each time step, in LISA frame 
(including Doppler shift correction to the frequency). The corresponding SNR at each 
time-step is referred to as instantaneous SNR i.e. snr+^x- We integrate the instantaneous 
SNR as given below 

SNR2^(ivb) = r sml^{eL{t;wB),(l)L{t;wB))dt. (9) 
Jo 

The network SNR while tracking wb over 1 year period is obtained by summing the 
squared integrated SNR's of the individual data streams ^ 

3 

SNRgfetCi^-B) = SNR2 (wb) + SNR2 (wb) = £ SNR?(ivb) . (10) 

1=1 

In Fig. 2(b), we plot the LISA sensitivity S = 5/SNRint for a monochromatic source 
tracked with several data combinations. SNRint is the integrated SNR along the source 
track (in LISA frame) for a given combination. For switching combination, e.g. optimal 
streams v+,x, SNRint = SNR+,x- 

Fig.2(b) displays the following features: (i) At low frequencies (/ < 3 mHz), v+, Vx 
have similar sensitivities and are proportional to (similar to y*^^'^^); the data stream 
(— ^) is insensitive to GW. (ii) Above 25 mHz, all F^^^'s become comparable in 
their sensitivities. The wavelength of GW (Agw) is comparable to L. This introduces 
geometry dependent features in the LISA sensitivity curve. For example, Y^^^ combina- 
tion is most sensitive when /g^ = fL = 1/^ = 60 mHz or multiples of f^. As V-|- contains 
contribution from Y^^\ this feature also appears in its sensitivity, see Fig. 3(a). The dips 
in the integrated SNR of v+ correspond to the zeros of y(3) (at / = (2n - 1)/l/2 mHz, 
n is a positive integer). 

3.2. Beaming of Optimal Stream 

Although, tracking wb with a network comprising of Y^^^^ — Nety — and with a 
network of v+ and Vx — Net+.x — gives the same integrated (as well as instantaneous) 
SNR [see, for instance, Eq.(lO) and Fig.2(b)], we show in this section that they possess 
completely different beaming properties. 

To demonstrate this, we divide the sky in the barycentric frame in 1 square degree 
patches. We observe each of this patch with the data streams 7*^^^ and integrate the 
squared SNR for a year along its track in the LISA frame (Hq = 1). In order to perform 
the same exercise with v+, Vx, vq, we choose to track an arbitrary source direction 
{Ob,^b} = {60°, 200°}. In Fig. 4,5(a), contour plots of integrated squared SNR are 
drawn (in barycentric frame) for all the data streams. The SNR is normalized with 
respect to SNR+. The observations and the implications are as follows: 



^ Note that SNR^g((M>B) can also be obtained by summing the squared SNR's of the individual data 
streams (trace of p^^jj) and integrating the resultant. 
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FIGURE 3. (a) SNR of v+ and Y^''> vs f, (b) SNR + /SNRNetVs/ for tracking {0b, = {60°, 200°}. 



(i) The data streams v+, Vx and vo are all beamed. In other words, the beam pattern of 
these streams peaks towards the tracking direction, whereas the Y^'^ streams are sensitive 
to a large fraction of the sky; see, for instance. Fig. 4(b), 5(a). The contours of the 
integrated SNR while tracking a particular direction give the point spread function (psf) 
of the source. In the language of the GW data analysis literature, these contours can be 
related to the ambiguity function in the source location parameter space. The width of 
the psf determines the size of the template. 

(ii) Net_|_ X gives the same integrated SNR as that of Nety. However, Net+^x is highly 
beamed as opposed to Nety. The complimentary feature of Net+.x and vq is apparent 
from the bright and dark patch centered around the tracking direction [see for instance. 
Fig. 4(b), 5(a)]. This property can have a possible immediate application in LISA data 
analysis. In LISA, we expect to observe a large number of GW sources from different 
sky locations in nearby frequency bins. Thus, with the combination of Net+^x and its 
complementary null stream vq, one could systematically suppress other sky directions 
without compromising on SNR. 

(iii) In the low frequency (/ < 10 mHz), Net+ x does not show beaming and gives the 
same beam pattern as that of Nety. This is because the angular response of y(i)'(2) is 
that of a single Michelson interferometer. The stream F^^) differs from the Y^^^ by 45° 
rotation [see [2],Fig.4 (a)] which makes the pattern azimuthal invariant after tracking. 
Further, the angular response of any TDI combination will be limited by the signal 
antenna angular pattern (size of earth's orbit ~ Agw)- 

(iv) In Fig. 4(b) and 5(a), v+ and Net+^x are beamed along the tracking direction 
{60°, 200°}. At 50 mHz (being closer to / ~ /l) , as opposed to at 25 mHz, y(3) 
more sensitive compared to 7'^^^''^^). 

One would expect that beam-width of the optimal streams is a monotonically decreas- 
ing function of frequency. However, we find super-imposed oscillations on the monoton- 
ically decreasing behavior in the beam-width[see. Figs. 5(b),6]'^. Below, we summarise 



In order to test this, we choose v+ (maximum SNR) to track the source at different frequencies. This is 
because, it is shown in Fig. 3(b),, for / > 25 mHz, v+ contributes more than 80% to Net+_x ■ 
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FIGURE 4. Integrated squared SNR at (a) / = 10 mHz, (b) / = 25 mHz and tracking {0b,<|)b} = 
{60°, 200°} with various data streams. 
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FIGURE 5. Tracking (^>ij} = {60°, 200°}: Integrated squared SNR (a) at / = 50 mHz for various 
data streams, (b)for v+ at different frequencies. 



the features of Fig. 6: At low frequencies (< 25 mHz), the overall beam-width of v+ 
decreases monotonically and is dominated by F^^'^^. At frequencies above 25 mHz, one 
observes modulation pattern in the beam-width (with frequency //,). More importantly, 
contributes a large fraction to SNR+ as compared to the other two F's, 
see Fig. 3(a). Hence, the beam-width of v+ is expected to carry the features of Y^-^\ 

As we know, Y^^^ is completely symmetric between the 3 arms. Its angular response is 
insensitive in the neighborhood of LISA's zenith as well as its equatorial plane whereas 
it is sensitive in the polar window of 6l '■ {30°, 60°}. A source at 6b = 90° follows 



a circular track in LISA sky with 6l = 60°, see Fig. 2(a). Hence, the source's track 
coincides with the sensitive part of 7^^^ beam-pattern which explains its maximum 
integrated sensitivity along Gg = 90°. As frequency increases, especially near / = nfi, 
due to the symmetry of Y^^\ for all Oi, the integrated antenna pattern of Y^^^ becomes 
invariant to azimuth 0b, see Fig.4(b),5(a) which results in increase in the beam-width. 
This explains the sudden increase of solid angle at nfi in Fig. 6. 

From Fig. 6, one can estimate the number of non-overlapping patches (templates) 
required to cover the entire sky. For instance, at / ~ 10 mHz, one requires 20 such 
patches whereas at / = 30 mHz, this number goes to 60. 
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FIGURE 6. Tracking {0b, 0b} = {60°, 200°} with v+: beam-width (obtained at 90% level) vs f. 



4. CONCLUSIONS 

In this work, we have studied the beaming properties of directionally optimal data 
streams v+,Vx and vq which contain the complete information of the GW signal in a 
particular direction. We have shown that they are beamed (after tracking the source for a 
year) and could be useful in eliminating certain sky directions or for consistency checks 
in LISA data analysis. 
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